Low temperature physics of Sr2Ru04: multi-component superconductivity 
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It is shown that the properties of the superconductivity 
observed in Sr2Ru04 can be accounted for by properly tak- 
ing into account the internal degrees of freedom of the system. 
The electrons are strongly correlated in the Ru02 planes while 
the superconducting phase transition is driven by "umklapp 
scattering" between the planes. The in-plane magnetic corre- 
lations are determined via a non-perturbative approach. The 
orbital degeneracy and the magnetic moments are internal 
degrees of freedom of the Cooper pairs that account for the 
thermodynamic properties in the superconducting phase. 



With the discovery of the high temperature supercon- 
ductors a whole class of transition metal oxides became 
a focal point in condensed matter research. These ma- 
terials exhibit many unconventional properties most of 
which are interpreted controversially to date. A number 
of experimental probes [|]-]3) show that the paired elec- 
trons in Sr2Ru04 carry a magnetic moment making its 
superconductivity unconventional 0| . The present paper 
shows that Sr2Ru04 is a striking example of how uncon- 
ventional physics can have conventional origins. 

Sr2Ru04 is the first layered transition metal oxide that 
exhibits superconductivity in the absence of copper ions 
1^ . The lattice symmetry is tetragonal and isostructural 
to La2Cu04 with lattice parameters a = b — 3.87 A in 
the Ru02 plane and c = 12.74 A out-of-planc. In con- 
trast to La2Cu04 the ruthenate is a Fermi liquid in the 
temperature range of Tc < T < 30 K [^] . Superconduc- 
tivity appears below Tc ^ 1 K. 

The significant correlations ||^,^ in Sr2Ru04, the S — 1 
moments on Ru'*+ impurities in Sr2lr04, and ferromag- 
netic correlations in SrRuOa led Rice and Sigrist to 
propose that the superconducting order parameter has 
p-wave symmetry promoted by ferromagnetic correla- 
tions analogous to ^He. A similar proposal was made by 
Baskaran [|lO| based on a comparison with high Tc ma- 
terials and emphasizing the role of Hund's rule coupling. 
This idea is supported by experiments that show that 
the static magnetic properties of Sr2Ru04 are the same 
in the normal and the superconducting phase while 
no conclusive experimental proof |pd]-|l3|| for the p-wave 
symmetry of the superconducting order parameter was 
found. No indication for ferromagnetic correlations has 
been found neither in neutron scattering investigations 
nor other approaches p5|-p7|. Furthermore, the spe- 
cific heat pS] and nuclear quadrupole resonance (NQR) 



[0 are consistent with two-dimensional gapless fluctua- 
tions in the superconducting phase of Sr2Ru04 which are 
absent in superfluid '^He. 

The present approach consists of three logical steps. 
First it is shown that the correlations in the RUO4 planes 
are dominated by quantum fluctuations. Based on non- 
perturbative methods p9[ | it is then possible to account 
for the experimentally observed magnetic structure factor 
^ in detail. 

In the second step it is shown that the coupling be- 
tween the RUO4 planes enhances the electron pair corre- 
lations via umklapp scattering processes ]2C| ] . In contrast 
to the in-plane correlations the inter-plane coupling can 
be treated perturbatively and leads to the superconduct- 
ing phase transition. The presence of the two degenerate 
and dyz orbitals together with Hund's rule coupling 
leads to Cooper pairs that carry a magnetic moment. 

In a third step it is shown that the approach is consis- 
tent with a multitude of experimental observations 
The magnetic moment of the Cooper pairs accounts for 
the static magnetic properties 0,1 of Sr2Ru04. Dis- 
cussed closer herein is the observed specific heat. 

Incommensurate magnetism In order to develop an 
understanding of the electronic correlations in the RUO2 
planes it is instructive to first discuss how perturbative 
approaches fail to describe the results from neutron scat- 
tering [ p^ in the normal phase. 

The bare electronic band structure of Sr2Ru04 has 
been determined unambiguously by a number of meth- 
ods Based on a tight binding band model for 
the band structure the effects of the on-site Coulomb re- 
pulsion can be included by a perturbative approach, the 
random phase approximation (RPA). It is then possible 
to calculate the magnetic structure factor. Various such 
studies have been performed with similar results [^2|-^ . 

Figure |l| visualizes a typical result for the imag- 
inary part of the dynamical magnetic susceptibility in 
RPA based on the band structures from the literature. 
The data have been convoluted with the experimental 
resolution appropriate for inelastic neutron scattering 
experiments. The structures in Fig. [l] reflect the dis- 
persion of the bands in the qx-qy plane in reciprocal 
space. While the results from neutron scattering show 
pronounced peaks at [±0.6^, ±0.6^] many of the other 
prominent structures have not been observed [ p^j26| . 

An obvious explanation is that the perturbative results 
underestimate the quantum fluctuations of the quasi two- 
dimensional system. In the presence of quantum fluctu- 
ations the magnetic susceptibility is expected to be en- 
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FIG. 1. Typical plot of the magnetic susceptibility in the 
qx-qy plane from RPA. Data are folded with the experimental 
resolution in neutron scattering experiments. 



hanced much more homogeneously. The question then 
arises about the origin of the experimentaUy observed 
peaks at [±0.6f , ±0.6f ]. 

The answer can be given in terms of the correlations 
of the dzx and dyz orbitals. The kinetic energy of the 
electrons in these orbitals is dominantly one dimensional 
p7| , p^ p|] , namely along x for the dzx electrons and along 
y for the dyz electrons. This allows us to bosonize the 
kinetic energy of the electrons in these orbitals. The 
on-site Coulomb interaction can then be included in a 
non-perturbative way. 

The Fermi velocity of the dzx and cL^ electrons has 
been estimated to be vp ~ 0.7 eVa jfJJT^. Perturba- 
tive approaches j2^j2^,|2^,Q involve values for the on-site 
Coulomb repulsion that are of the same order as v-p, i.e., 
the system is rather in an intermediate than weak cou- 
pling regime. The coupling hybridizes the dzx and dyz 
electrons leading to dominantly quasi one-dimensional 
magnetic correlations along the diagonals of the system, 
i.e., along x = ±y or qx — The mathematical 

derivation of those correlations from a microscopic model 
is rather involved and the interested reader is referred to 
a forthcoming article (Ref. JTot). 

The resulting anticipated quasi one-dimensional el- 
ementary magnetic excitation spectrum in reciprocal 
space is depicted in Fig. |[ The spectrum is gapless at 
(±[^-2fcF],±[^-2fcF]) = (±0.6f,±0.6f) which leads 
to a dynamical magnetic structure factor that is strongly 
peaked as observed experimentally |Q . kp is the Fermi 
wave number of the dzx and dyz electrons. The existence 
of a continuum is also consistent with experiments. Most 
striking though is that the width of magnetic excitations 
at energy lo 
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is basically independent of the temperature as observed 
experimentally iQ (red bar Fig. H). This feature is diffi- 
cult to explain within other theoretical approaches. The 
comparison of Eq. (0) with experiment yields a value of 




FIG. 2. Effective, quasi one-dimensional magnetic excita- 
tion spectrum in the qx-qy plane as anticipated from the dzx 
and dyz electronic subsystem at intermediate coupling. 



the velocity of the elementary magnetic excitations of 
« 120 K which is consistent with the temperature 
dependence of the peak intensities ||l^,|2^ • 

In conclusion the intermediate coupling electronic cor- 
relations in the Ru02 planes of Sr2Ru04 must be de- 
scribed by a homogeneous, two-dimensional part and 
quasi one-dimensional magnetic correlations along the di- 
agonals of the basal plane. 

Pair instability For the system to undergo a phase 
transition the correlations in the planes have to be cou- 
pled to form a three dimensional system. The inter-plane 
hopping t± between the dzx and dyz orbitals has been 
estimated to be an order of magnitude smaller than the 
in-plane hopping [p7|J^ . We thus have a small parameter 
justifying a perturbative approach. The inter-plane hop- 
ping involving dxy orbitals is expected to be significantly 
smaller for geometrical reasons [p|j27|,^ and is neglected. 

Let us consider the lowest order coherent inter-plane 
pair hopping term 
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The electron creation and annihilation operators are 
^\ V a ^^"^ u a f*^^ orbital v S {x, y} with spin a E {1,1} 
on site I. The orbital indices v^^^iy^fi' € {x,y} are re- 
stricted by the Pauli principle for a' — a to fi' ^ fi and 

I/' ^ V. 

The prefactor 1 > ff^^^/ > takes into account 

that the pair hopping into the same orbital is suppressed 
by the Pauli exclusion principle. The presence of Hund's 
rule coupling favors spin triplet pairs over spin singlet 
pairs. Consequently one expects mixed orbital, spin- 
triplet pair hopping to be enhanced with respect to mixed 
orbital, spin-singlet pair hopping and same orbital hop- 
ping 

The Fermi operators in Eq. (^) can be Fourier trans- 
formed via 



(3) 
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The real-space positions of the Ru ions are given by 
Ri, the fc are wave vectors. Limiting the discussion to 
the channel with the largest pair hopping amplitude the 
Hamiltonian becomes in reciprocal space 

1: PtMPt,,{<i) ■ (4) 

b—x.y,z q 

The mixed orbital, spin-triplet pair operators are defined 
by 

Pt,b{l) — ^ '^^a,a" ^ct",(t' ^q-k,x,tT^k,y,iT' j (5) 

where cr^ ^, are the Pauli matrices with b G {x, y, z}. The 
effective potential is given by 



t i a a c 

Vq,t = — cos -q^ cos -qy cos -q^ 

tip z z z 
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For q = q with 



f 27r ^ 2t: , 27r „ 27r . 27r . 27r J 

qe<^ — yH z, — x, — y, — z\ (7) 

y a a c a a c ) 

the potential is extremal and attractive, i.e., lq.,t 
—t\/v-p- As a consequence of the body centered tetrag- 
onal lattice structure the are reciprocal lattice vec- 
tors. In other words the body centered tetragonal lattice 
symmetry allows for an umklapp scattering driven pair 
instability. 

The phase transition occurs at T = Tc when the Stoner 
criterion 
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is fulfilled. Xtb is the in-plane pair-pair correlation 

function with the pair operators as given in Eq. (I). The 
strength of the effective potential can be estimated |2^] 
as —Vq.^t = Vo,t ^ t\/v-p « 6 K which is consistent with 
the observed small critical temperature Tc ~ 1.5 K. In 

the absence of a magnetic field and neglecting symmetry 

(p\ (p\ (p\ 
breaking effects Xt x — Xty — Xt z ^^'^ below the phase 
transition all three order parameter components have a 
finite expectation value (Pt,{,(0)) ^ V &. 

The superconducting phase transition is induced by 
inter-plane umklapp scattering. The Cooper pairs carry 
spin one that is the combined magnetic moment of the 
two bound electrons. This is consistent with the absence 
of a change of the Knight shift ^ and of the magnetic 
susceptibility ^ in the superconducting phase. 

Order parameter fluctuations The order parameter 
in Sr2Ru04 has multiple internal degrees of freedom. 
The magnetic degrees of freedom are given by fluctu- 
ations between the three order parameter components 
{Pt.b=x.y.z) ■ They can be parameterized by a SO (3) vec- 
tor as Jig = (sin Os sin 6z , cos 6^ sin 9^ , cos . 
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FIG. 3. Specific heat. The dash-dotted line is a fit to the 
region of quadratic temperature dependence. The dashed line 
is the linear fit from Eq. (^^. 

Muon spin relaxation (/iSR) measurements and the 
out-of-plane critical field suggests that the magnetic 
moment of the Cooper pairs lies predominantly in the 
x-y plane The vector Vis is then reduced to S0(2) 
symmetry, i.e., 9z — 7r/2. This does not alter the physical 
results discussed herein. 

The mean-field approach as discussed in the previous 
section ignored possible spatial anisotropics of the order 
parameter. Critical field measurements are consistent 
with the presence of two spatial order parameter com- 
ponents [ ^J2l| . The model discussed in the context of 
the incommensurate magnetic fluctuations suggests that 
there arc two distinguished directions in the system. As 
becomes apparent in Fig. |^ those are along q^ — qy and 
qx = —qy. Indeed, it follows from the microscopic model 
that the order parameter has two spatial components as- 
sociated with these directions |^ . The components can 
be parameterized by a two-component vector S7f . 

In the superconducting (ordered) phase the fluctua- 
tions in each of the channels /i G {s,f} can be described 
by a non-linear sigma model with the action [pl[ 



drj: 



(9) 



The correlations in Sr2Ru04 along the z axis are two 
orders of magnitude smaller than in the Ru02 planes 
1^,111 . The stiffness p can thus be assumed to be isotropic 
in the plane but negligible along z. The action Eq. ^ is 
quasi two dimensional. 

At low temperatures the action Eq. (^) describes quasi 
two-dimensional, gapless fluctuations. Their specific heat 
is quadratic in temperature. The symbols in Fig. ^ 
show the experimental specific heat over temperature 
C/T as published in Ref. together with a linear low- 
temperature fit (dash-dotted line). 

For not too small magnetic fields the magnetic fluc- 
tuations described by Sg in Eq. (^ are suppressed and 
will not contribute to the specific heat. One expects the 
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slope of Ch/T versus T to be smaller than in the case 
without magnetic field consistent with the experimental 
observation 

A less obvious consequence of the gapless fluctuations 
in the superconducting phase is the linear temperature 
dependence of the specific heat below the phase transi- 
tion. The blue symbols in Fig. || show the data from 
Ref. [|l8| in the linear range near Tc together with a lin- 
ear fit (dashed line). The origin of this linearity can be 
understood by expressing the free energy of the system 
in terms of the electronic excitation gap A ^ \{P)\ in 
the superconducting phase Q valid for small reduced 
temperatures t = 1 — T/Tc ^ 1: 

Fa = Fo - A < + £) + S + 0{lS.^) . (10) 

The presence of the third order term can be motivated by 
integrating out the internal degrees of freedom fig and ilf 
|pO| . Since Eq. (^^ is an expansion in A > and not in 
the order parameter itself the third order term leads to a 
phase transition of third order in the sense of Ehrenfest's 
definition. 

The temperature dependence of the gap is obtained by 
minimizing the free energy Eq. ([L0[) as 



(11) 
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The linear temperature dependence of the gap at the 
phase transition has been observed via the ^SR rat e [ppO| 
and the excess current in tunneling experiments pl|20|. 
In a next step the specific heat near the phase transition 
C/T = -(92Fa)/(5A2) is found as 



C's|t<l _ Cn 
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(12) 



that is consistent with the linear temperature dependence 
at the phase transition. C„ is the specific heat in the 
normal phase. 

In conclusion it has been shown that the in-plane corre- 
lations of the Ru02 subsystem of Sr2Ru04 exhibit domi- 
nant quantum fluctuations making the application of per- 
turbative approaches questionable. The magnetic corre- 
lations have quasi one-dimensional components along the 
diagonals of the basal plane of the crystal. The super- 
conducting instability is induced by inter-plane umklapp 
scattering in the body centered crystal structure. The 
hopping amplitude is largest in the mixed orbital, spin- 
triplet channel. Possible competing long range magnetic 
order is suppressed because of the strong in-plane quan- 
tum fluctuations and the frustrating body centered crys- 
tal structure. 

The number of experimental results that are consistent 
with the model for the superconductivity in Sr2Ru04 pre- 
sented here goes beyond what can be discussed in this 
overview. The interested reader is referred to Refs. ||2C| ] 
and pl|. The new insight into the superconductivity in 



Sr2Ru04 may well help clarify the open questions posed 
by the properties of other layered transition metal oxides. 
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